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IT) ' Abstract 

00 ■ 

■t;:j|- ■ We apply the resummed version of the Liischer formula to ana- 

lyze finite volume corrections to the mass of the nucleon and of heavy 
•/^ ' mesons. We show that by applying the subthreshold expansion of the 

^^ ■ scattering amplitudes one can express the finite volume corrections 

in terms of only a few physical observables and the size of the box. 
In the case of the nucleon, the available information about the quark 
mass dependence of these physical quantities is discussed and used to 
assess the finite volume corrections to the nucleon mass as a function 
5^ I of the quark mass including a detailed analysis of the remaining un- 

certainties. For heavy mesons, the Liischer formula is derived both 
fully relativistically and in a nonrelativistic approximation and a first 
attempt at a numerical analysis is made. 
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1 Introduction 

Lattice calculations have recently almost reached physical pion masses (see 
e.g. Refs. [1, 2]), thereby establishing firm contact with chiral perturbation 
theory (ChPT) [3-5]. The effective theory therefore becomes increasingly 
useful to address systematic effects in lattice calculations like the finite vol- 
ume, for example. To use the effective chiral Lagrangian to evaluate such 
effects was first advocated in Refs. [6, 7] where it was also shown that in 
the p-regime, where the Compton wavelength of the pion is much smaller 
than the size of the box, the infinite volume Lagrangian can be used for 
doing calculations. An alternative approach to evaluate finite volume cor- 
rections (FVC) to masses in the p-regime is the Liischer formula [8, 9]. This 
relates the FVC of a particle P to an integral over the forward scattering 
amplitude of the same particle P off the lightest particle in the spectrum. 
The main contributions to this integral come from the low-energy region, 
while high-energy contributions are exponentially suppressed. This justifies 
the use of amplitudes evaluated in ChPT as an input. An improved version 
of the Liischer formula, where subleading terms are resummed, was already 
successfully apphed to meson masses and, in a generalized form, to coupling 
constants [10-13]. FVC for the nucleon mass in ChPT and the Liischer 
formula have been discussed in [14-16]. 

One of the advantages of the Liischer formula is that it allows a relatively 
easy estimate of higher order effects and yields simple analytical expressions. 
In Ref. [12] it was shown that inserting the tree-level tttt scattering amplitude 
into the resummed Liischer formula exactly reproduces the one-loop result 
of the FVC for the pion mass evaluated in ChPT. At the two-loop level, this 
does not hold anymore, because in the derivation of the Liischer formula, only 
one propagator is considered to be in finite volume — two-loop contributions 
in which propagators appearing in different loops are simultaneously taken 
in finite volume are exponentially subdominant and neglected from the start 
in Liischer's approach. However, the two-loop FVC of the pion mass were 
explicitly calculated and compared with the result of the resummed Liischer 
formula with a one-loop amplitude in the integrand in Ref. [17]. It was found 
that the difference is very small and thereby a solid argument was provided 
to base further analyses of FVC on the resummed Liischer formula. 

In this paper we apply the latter to estimate the FVC to masses of heavy 
particles, and, in particular, to nucleons and to mesons containing a heavy 
quark. Although these two kinds of particles are very different — nucleons are 
spin-y2 particles and are much lighter than the spin-0 B mesons, for example 
— the essential common feature is that both are much heavier than the pions. 
Neither the exact value of their mass, nor their spin, make a qualitative 
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Figure 1: Singularities of the forward scattering amplitude in the complex 
i^-plane. The cut on the real axis is due to the intermediate state Hn and 
the pole on the left-hand side of the imaginary axis to the intermediate H* 
state (or H in the case of the nucleon). The position of the pole depends 
on m,'jj, — m,'jj and moves to the right-hand side of the imaginary axis if 
m,^fj* — m,'jj > M^. Notice that the amplitude is an even function of u and 
that for clarity we have drawn only half of the singularities. 



difference in the evaluation of the FVC What makes the treatment of the 
two cases very similar is that the singularity structure of the ttH scattering 
amplitude at low energy [H being the heavy particle) is dominated by a single 
pole due to the exchange of a particle degenerate (or almost degenerate) with 
H. The forward scattering amplitude is a function of only one kinematical 
variable u. The pole due to the exchange of a particle H* in the s-channel is 
located at 



z/ = -(M^ - ml. + m]j)/{2mH] 



mjjt^mH 



-Mll{2mu) , 



(1) 



as illustrated in Fig. 1. If m'jj, — rn^ > M^ the pole is on the right-hand 
side of the integration path (which coincides with the imaginary axis) and 
the Liischer formula contains only the term with the integral. As we take the 
limit ruH* — ;■ Tan the pole moves to the left until it crosses the imaginary axis 
and finally reaches u = —M^/{2m,H)- When the pole crosses the imaginary 
axis an extra local contribution to the Liischer formula appears [8], which 
is exponentially leading. The study of the heavy meson case, in which we 



can vary A^, = rriH* — rriH and even send it to zero, allows us to better 
understand the role of the local, exponentially leading contribution in the 
Liischer formula. We stress that this singularity structure is different from 
the one in the case of the light pseudoscalar mesons L considered in Ref. [12]. 
In that case the most important low-energy singularity is the cut due to 
the 7tL intermediate state (also shown in Fig. 1 starting aX u = M^^ for all 
particles): the closest single pole is due to the exchange of a vector meson, 
but is at higher v values. 

For the case of the nucleon we can rely on existing calculations of the 
infinite volume vriV scattering amplitude to one loop and on the knowledge 
of the low-energy constants (LEG) appearing therein. In particular we show 
how the application of the so-called subthreshold expansion — very well 
known in the phenomenology — leads to a particularly simple formula for 
FVC of the nucleon mass. The solid phenomenological knowledge of the 
subthreshold coefficients allows a very precise determination of the FVC at 
or near the physical value of the quark mass. In our numerical analysis 
we show how the error in the evaluation of the FVC increases as we go to 
higher quark masses and identify the LEC which control these effects. In 
particular a better understanding of the quark mass dependence of the axial- 
vector coupling constant of the nucleon would lead to a much better control 
of the FVC. The relevance of this issue for the calculation of the FVC for 
the nucleon mass has been already discussed at length in Ref. [K)], where we 
presented preliminary results (dating already five years back). Indeed until 
lattice calculations of qa at lower quark masses were published, we were not 
able to present a reliable numerical analysis of these FVC. The situation 
has now clearly improved, but is not yet fully satisfactory and we do think 
that there is still a lot to be understood about the behavior of qa and other 
nucleon-related quantities as functions of the quark mass. 

Only very few publications on finite volume effects for heavy mesons exist 
so far [18-22], most of which are not concerned with the FVC to the mass. 
Also about the scattering amplitude of pions off heavy mesons a lot less is 
known, even only in the low-energy region. We give the expression of the 
tree-level scattering amplitude and show how the Liischer formula needs to 
be modified if one considers the heavy mesons away from the infinite mass 
limit {i.e. as nondegenerate) . We show that also in this case one can define a 
subthreshold expansion completely analogous to the case of the nucleon, and 
express the FVC in terms of the subthreshold parameters, so emphasizing the 
close analogy between the two cases. We provide a simple numerical analysis 
on the basis of the limited information which is available on the irHH* (for 
H = D or B mesons) coupling constant and discuss the difference between 
the fully relativistic case and the nonrelativistic limit of the Liischer formula. 



2 Resummed Liischer formula for the nucleon 
mass 

In this section we present and discuss the formula we will use for analyzing 
the finite volume effects on the nucleon mass numerically. We rely on the 
resummed version of the Liischer formula [8, 9] which has been proposed 
in Ref. [12]. The formula for the relative finite volume correction i?7v = 

(mAr(L) — mAr)/mAr reads as follows: 
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where 



e^ = -^, K = M^L, D+{y) = m^D+{iM^y,0) (3) 

ZrriN 

where D'^ is one of the components of the elastic vrA^ scattering amplitude, 
which is defined as 

T{7r''iq)N{p) ^ n'^' iq')N{p')) = T^^a = Sa'aT^ + l[ra',r,]T- . (4) 

Each of the two isospin components is then broken down into 

(5) 
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and each of the amplitudes depends on the kinematical variables t and u, 
defined as 

t = (q — q'Y , u = , where s = (p + qY , u = (p — q'V . (6) 

4m7v 

The term with n = 1 in Eq. (2) has been first given by Liischer [8] as an 
application of his general formula for finite size corrections to masses [9], but 
with a factor 2 missing in the first term within the brackets. This has been 
pointed out in Ref. [11] and later confirmed in Ref. [15] — it is indeed easy 
to follow the steps in Liischer's general proof [9] adapting it to the present 
case and derive formula (2) as it stands. 
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Table 1: The multiplicities m^n) for n < 15. 



The extension to the exponentially suppressed terms with n > 1 can 
be immediately obtained following the general derivation given in Ref. [12]. 
They appear weighted by the multiplicities 771(71) — these are easily calculable 
and are reproduced here for convenience in Table 1 for n < 15. 

The resummed version of the formula exactly reproduces the one-loop 
result in ChPT if the tree-level scattering amplitude is inserted in the integral 
in (2), see Ref. [12]. In the present case, due to the special counting rules 
in the one-nucleon sector, the ttN amplitude receives tree-level contributions 
both at order p as well as at order p^. Inserting the irN amplitude up to order 
p^ one can indeed exactly reproduce the calculation of the 0{p^) and O(p^) 
finite volume corrections performed in Ref. [11]. Since the ttN amplitude 
is known well beyond the 0{p'^) level, in fact up to O(p^) [2:>], one can use 
this formula to go beyond the one-loop level. In this manner the corrections 
beyond 0{p'^) will be given exactly as far as the first three leading exponential 
terms are concerned, and only approximately for the subleading exponential 
terms beyond e~^^^'"^ [12]. Based on the analysis in Ref. [17] we expect the 
contributions which go beyond the resummed formula (2) to be numerically 
small. 

As seen in (2,3), the irN amplitude is needed here in a particular kine- 
matical configuration, namely for t = and for v purely imaginary and small: 
the contributions with large values of v are suppressed by the exponential 
kernel in the integral in (2). It is therefore natural to make a Taylor expan- 
sion of the amplitude around z/ = after having subtracted the pole due 
to the one-nucleon exchange diagram (also called the Born term). Such an 
expansion is in fact already well known in the phenomenology and is referred 
to as the subthreshold expansion. It reads as follows: 

D+{u, 0) = Dl{u, 0) + D;{u, 0) + DtAs. u) (7) 

where 

^^v(-,0) = ^V^ 
^ rriN ^B — ^ 

d;{^.^) = 4o + rf>' + 4o^' (8) 

and ub = —Ml/{27njq). The function D^^{s,u) contains the analytically 
nontrivial part of the amplitude. Up to order p'^ in the chiral expansion this 
can be written as a sum of two single variable functions: 

D+ (s, u) = Dfis) + Dtiu) + 0(/) (9) 

which admit the following dispersive representation: 

^1 s = — / dv'-—-^^^— 10 



where s' = 2??7,jvz/' + mj^ + M^. This representation shows that, due to the 
large number of subtractions, the function Df is small near z/ = 0. As we 
will see later, its contribution to the finite size shift of the nucleon mass is 
negligible. 

This observation leads to the following expression for the relative finite 
volume shift Rn: 
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where 



-y(l+y2)A. 



/pv(A.,£.)=/ dy , B\K) = / dy y'e~V(^^y'^>^^ (12) 

are the relevant finite volume integrals, J^ = nij^M^^d^Q and -RAf,na is the 
remainder coming from the (subtracted) analytically nontrivial part of the 
amplitude. 

If we neglect the contribution RN,na the representation (11) is very simple 
and expresses the finite volume shift of the nucleon mass in terms of only a 
handful of physical observables: the pion and proton masses, M^r, m^r, the 
pion-nucleon coupling constant g^^N and the three subthreshold parameters 
d^. If one knows the low-energy constants which appear in the chiral repre- 
sentation of these quantities, one can predict their quark mass dependence 
and therefore the finite volume shift Rn as a function of the quark mass. 
We also stress that for simulations performed at the physical point, Eq. (11) 
allows an evaluation of the FVC based only on input extracted from the 
phenomenology, and indeed on rather well-known quantities like Qt^n and the 
subthreshold coefficients. 



3 Quark mass dependence of ttiat, qa^ Q-kN and 
the subthreshold coefficients 

The chiral representation for the observables ttitv, Q-kN up to order p'^ can be 
found in Refs. [23, 24] and will be reproduced here for convenience. It reads 



as follows: 
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where A^ = IGtt^, the symbols m, g and F indicate the nucleon mass, the 
axial coupling constant and the pion decay constant in the SU(2) chiral limit, 
whereas M^ = 25m is the leading term in the expansion of the pion mass 
square in powers of the quark mass. The LEG Cj, di and Cj appear in the 
chiral Lagrangian in the one-nucleon sector at order p"^, p^ and p^ respectively, 
and the tildes indicate renormalized, scale independent versions thereof [23] 
(the Cj's as well as dig are scale independent — see Appendix A for more 
details). 

For the axial charge, the order p^ result^ is partially known [25] and will 
be used in the following. 
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Note that the above expressions for 74 and (3^ do not represent the full 
0{p^) result. We introduce a generic LEG / which collects all the polynomial 
contributions proportional to M^. The chiral expansion of the subthreshold 
coefficients reads: 



d+ 
"00 



2(2ci-C3)M2 ^ g^3 + 8g')Ml 

647rF4 

3 ig' + Qg') 



F2 



+Mt 



Cl 



es 



FJ Svr^Fi 



% + 



647r2Fi 



m 



2ci - C3 
167r2Fi 



^We refer here to the counting in the Lagrangian. As seen in Eq. (14) this translates 
into an expression including up to the 0(A/^) correction in gA- 
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3.1 Determination of the LEG 

The LEG from the nucleoli sector are determined by fitting the chiral rep- 
resentations of rrifyf, qai gnN and df^ to experimental data and results from 
lattice simulations (on rriN and qa only). Details of the fits are given in 
Appendix B. 

In order to have a better handle on the uncertainties, we perform a fit 
to the order p^ chiral representations as well as to the complete formulas 
indicated in Eqs. (13), (14) and (16). Note that in the following, we give 
the value for the scale independent Jig evaluated at the physical pion mass 
(see Appendix A) while for the higher order LEG we indicate the values of 



<(/^),/(/^) at /i = Mp. 






• At order p"^, we obtain 






m = 0.896 ± 0.003 GeV 


9 = 


= 1.126 ±0.007 


ci = -0.54 ± 0.04 GeV"^ 


C2 = 


= 1.79 ± 0.03 GeV"^ 


C3 = -3.37 ± 0.18 GeV"^ 


di6 = 


= 0.06 ± 0.11 GeV~2 


dis = -1.13 ± 0.19 GeV-2 







xVd.o.i. = 5.7/7. (17) 

Note that the term proportional to M^ in g^ is also included in the 
order p"^ fit, although it is a term of order p^. Furthermore, for the 
subthreshold coefficients, only the leading terms of the chiral expansion 
are retained at order p^. In particular, d^Q = 0. 

The fit with the untruncated expressions of Sec. 3 does not allow one to 
determine all of the appearing LEG. We have checked that if one holds 
three of the LEG at a constant value then the fit becomes stable. We 
choose to fix the three LEG C2, C3 and C4. To assess the uncertainties 
thereby introduced into the analysis, we use two different sets of values 



for these three LEG, given in Ref. [23], 

Setl C2 = 1.7GeV"^, cg = -3.6 GeV"^ , C4 = 2.1GeV"^, 
Set II C2 = 2.7GeV"^ , cg = -4.5 GeV"^ , C4 = 2.4GeV-^. 

In the following, this fit will be referred to as the order p^ fit, although 
the formulas also contain terms of higher order. For set I, the fit yields 



m = 


0.907 ± 0.021 GeV 


9 = 


- 1.201 ±0.077 




Cl = 


-0.56 ± 0.30 GeV^^ 


die = 


: -3.05 ± 1.54 GeV^ 


-2 


di8 = 


-1.20 ± 0.22 GeV-2 


f = 


:3.2±47GeV-^ 




e\ = 


:15±13GeV~=^ 


el = 


: -34 ± 72 GeV"3 




e^ = 


:37±16GeV~^ 


el = 


: 23.2 ± 1.4 GeV^^ 




xVd.o.f. = 


: 4.8/5 






(18) 


e for set II 


we find 








m = 


0.905 ±0.021 GeV 


9 = 


: 1.207 ±0.077 




Cl = 


-0.60 ± 0.30 GeV"^ 


die = 


: -3.43 ± 1.56 GeV^ 


-2 


di8 = 


:-1.2±0.2GeV-2 


f = 


:-0.1±47GeV-2 




e\ = 


:15±13GeV-3 


el = 


: 56 ± 72 GeV^^ 




el = 


: -86 ± 20 GeV"^ 


el = 


: 23.2 ± 1.4 GeV"^ 




xVd.o.f. = 


: 4.8/5. 






(19) 



The errors indicated for the LEG are statistical only and do not account, in 
particular, for the uncertainties due to higher order effects. The statistical 
error quoted above most likely underestimates the real uncertainties. 
Gomparing the different resulting values of LEG, one observes that the nu- 
cleon mass in the chiral limit lies very close to 0.9 GeV in all three fits. Also 
the values for ci end up quite close together, but fall out of the error bars 
for Cl quoted in Ref. [26]. However, given the fact that the Cj are correlated, 
we should compare with the value of Ci given in Ref. [23], which is in better 
agreement. Our fits with the two sets of order p^ LEG lead to consistent re- 
sults for Cl. The LEG of the order p^ chiral Lagrangian are less well known. 
Gomparing with Ref. [27], the values of c/ig of the order p^ fits agree remark- 
ably well. In the order p"^ fit, the resulting value of diQ is not reconcilable 
with phenomenology (see also Ref. [16]). We come to the same conclusion as 
the authors of Ref. [25]. Only the inclusion of the partial order p^ corrections 
in gA leads to a reasonable fit. The coupling dis is in agreement with the 

10 



range of values found in Ref. [28]. Since dig is determined with the help of 
the Goldberger-Treiman relation at the physical point, different values of dig 
merely reflect different values oi g. As for the LEG of the order p^ chiral La- 
grangian in the nucleon sector, the huge statistical errors of the fits confirm 
that they are basically unknown. 

3.2 Numerical analysis 

Given the values of the LEG from the fits, we plot the quark mass depen- 
dence of the quantities rriN, qa and the subthreshold coefficients. As can be 
seen in Fig. 2, the nucleon mass is described remarkably well by the order 
p^ curve which is just a quadratic function in the pion mass. Figure 2 also 
shows the axial charge. It exhibits the mild quark mass dependence observed 
in lattice calculations already at order p^ and stays within the order p'^ error 
bars over a wide range of quark masses. However, this picture hides the 
fact that the chiral series of qa is not well behaved, as the observed smooth 
behavior is the result of compensations among different chiral orders. We 
believe, however, that it is rather unlikely to have a strong quark mass de- 
pendence in the region between the physical pion mass and the present lowest 
lattice values (around M^ ~ 0.3 GeV). As we will show in the next section 
the present knowledge on the quark mass dependence of qa (which is based 
mainly on lattice calculations, rather than on chiral predictions) does allow 
us to calculate the FVG for the nucleon mass sufficiently reliably. In the 
latter it is actually the pion nucleon coupling constant which appears and 
Fig. 3 shows that its quark mass dependence is rather weak and that the 
difference between the two sets is negligible. 

The behavior of the subthreshold coefficients doo ^^^ '^lo beyond the phys- 
ical point is not predicted by chiral symmetry. Going from set I to set II leads 
to substantial changes in the higher order LEG and therefore to a very dif- 
ferent behavior as a function of the quark mass. However, the extrapolation 
towards the chiral limit seems to work quite well. The coefficient d^Q seems 
to be in much better shape as it does not depend on the Cj. Set I and set 
II almost yield identical results. However, this nice agreement should not 
mislead one into thinking that everything is under control here: large higher 
order corrections may modify this nice picture, since Eq. (16) is the leading 
order result for d^o- For instance, replacing F-n- by its value in the chiral limit 
in d^o yields the dashed-dotted curve in Fig. 4. Note that the coupling eg 
was readjusted to reproduce the correct physical value. 
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Figure 2: The upper panel shows the quark mass dependence of the nucleon 
mass and the lower panel of the axial charge, respectively. The dotted lines 
show the results of the order "j? fit, the solid lines represent order jy" set I 
and the dashed lines are order jr" set II. The error bands are represented by 
thinner lines plotted with the same line style as the central value. The circles 
show the lattice result as well as the physical point used as input. 
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Figure 3: The upper panel shows the quark mass dependence of the pion 
nucleon couphng constant and the lower panel of the scaled subthreshold 
coefficient M^^d^Q, respectively. The dotted fines show the results of the order 
p^ fit, the solid lines represent order p^ set I and the dashed lines are order 
p^ set II. The error bands are represented by thinner lines plotted with the 
same line style as the central value. The circles show the physical value used 
as input. 
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Figure 4: The upper panel shows the quark mass dependence of the scaled 
subthreshold coefficient M^dl^ and the lower panel of M^rf^Q, respectively. 
The dotted lines show the results of the order p^ fit, the solid lines represent 
order p^ set I and the dashed lines are order p^ set II. The error bands are 
represented by thinner lines plotted with the same line style as the central 
value. The circles show the physical value used as input. The dashed-dotted 
line shows the rescaled coefficient M^d^g with F^, replaced by F, as described 
in the main text. 
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4 Numerical analysis of the FVC to the nu- 
cleoli mass 

With the quark mass dependence of all the quantities appearing in Eq. (11) 
established, the finite volume effects can now be calculated for pion masses 
below ~ 0.4 GeV and box sizes L > 2 fm. In Fig. 5 we show the finite 
volume effects for boxes of spatial extent L = 2, 3 fm. To illustrate the size of 
the different contributions to the FVC, we break down the corrections into 
the different terms of Eq. (11). This is shown in Fig. 6. The contributions 
from the Born term dominate the finite volume effects and are practically the 
same for all three fits. Since the quark mass dependence of the subthreshold 
coefficients is almost unknown, their contributions to the finite volume effects 
vary a lot. The shift from the order p"^ fit to the order p^ fits is mainly due to 
a vanishing contribution from (i^Q at order p^. Removing this shift by setting 
^20 to its physical value brings the order p"^ curve down into the error bars 
of the order p^ fits. 

At order p^, the fit with set I shows cancellations between d^Q and d^g- How- 
ever, these cancellations are absent for set II, where cIqq and df^ contribute 
very little but the contribution from d^Q stays the same as in set I. This leads 
to the large difference between set I and set II. The central values shift by 
twice the error bars. It is obvious that a precise prediction of the finite size 
effects is very difficult away from the physical point. One also observes that 
the finite size effects strongly depend on the values of the LEG q, i = 2, 3, 4 
through the subthreshold coefficients. Therefore, if the finite volume effects 
are well determined by lattice simulations, this would provide useful con- 
straints and possibly a determination of (some combinations of) the q. 

If we expand our formula for Rjsf and keep in the scattering amplitude 
only the 0{p'^) contributions we agree exactly with the one in Ref. [11]. A 
comparison of the numerics at this order is therefore uninteresting (since 
by inserting the same values for the LEG one trivially gets the same). We 
stress, however, that as our plots in Fig. 5 indicate, by truncating the chiral 
expansion at NLO one may get the misleading impression of having rather 
moderate uncertainties which do not grow very fast with the pion mass. 
The application of the resummed Liischer formula allows one to go to higher 
orders and to verify explicitly that: 

• corrections beyond NLO are important; 

• the uncertainties grow rather quickly with the pion mass, so that going 
beyond pion masses of about 0.3 — 0.35 GeV they become of order 

100%. 
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Table 2: Our final results for Rn. In the first column, the value of the pion 
mass is given and in the subsequent columns, the values of Rn for three 
different sizes of the box are listed. The star indicates that for these values, 
M^L < 2. 



We refrain from comparing our error analysis with Ref. [29], where FVC 
including error bars were calculated at one loop in ChPT. In Ref. [29], plots 
of the FVC are only given for values of pion masses and box sizes in a region 
which we consider to lie outside the region of validity of the effective theory. 
Our final results are shown in Table 2 where we give Rn for different 
values of pion masses and box sizes. The quoted central value is the average 
of the central values of set I and set II and the error bars are chosen such 
that they cover the error bands of both set I and set II. 



4.1 Physical point 

As we have discussed in detail in the previous sections, the main difficulties 
in the numerical evaluation of the FVC for the nucleon mass arise from the 
limited knowledge of the quark mass dependence of the observables which 
appear in the resummed Liischer formula. Since several lattice collaborations 
are reaching or plan to reach the physical value of the pion mass, it is inter- 
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Figure 5: The FVC Rn as a function of the quark mass for L = 2fm (upper 
panel) and L = 3fni (lower panel). The dotted lines indicate the result of 
the order p^ fit, the solid lines are the order p^ fit with set I and the dashed 
lines the order p^ fit with set II, respectively. Note that for L = 2 fm and 
L = Sim, the pion mass should fulfill M^ > 0.2 GeV and M^ > 0.13 GeV, 
respectively, in order to stay in the region of validity of the Liischer formula. 
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Figure 6: Individual contributions to the Liischer formula (shown by lines of 
different colors, see legend). The finite volume effect Rn (black lines) is the 
sum of all the contributions. The line style indicates the order and the set 
of LEG used. The dotted lines denote the result of the order p^ fit, the solid 
lines are the order p^ fit with set I and the dashed lines the order p^ fit with 
set II, respectively. The size of the box is L = 3 fm and therefore, the pion 
mass should fulfill M^^ > 0.13 GeV in order to stay in the region of validity 
of the Liischer formula. 

esting to evaluate the FVG for this case — the advantage of our approach is 
that at the physical pion mass we can insert directly the measured values for 
the different observables. 

In Fig. 7, Rn is plotted as a function of the size of the box L. The 
plot shows the order p^ curves as well as both order p^ curves. A difference 
is visible in Rn only because at order p^, the subthreshold coefficient d^g 
vanishes. The main contribution to the FVG is given by the Born amplitude. 
The subthreshold coefficients all yield corrections of similar size. In fact, 
the contributions from df^ and d^g cancel to a large extent. Since the error 
bands are constrained very much by the small errors of the experimental input 
parameters we can predict the FVG very accurately for lattice simulations 
with physical quark masses in a box of the size L > 3 fm. If nucleon mass 
calculations with physical quark mass are done in a box of the size L > 3.5 fm 
one does not have to worry about finite volume corrections. The effects stay 
below 1 %. 
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Figure 7: FVC for physical pion mass as a function of the size of the box. 
The black lines denote Rn- The dotted line is the 0{p^) fit, the solid line and 
the dashed line are set I and set II of the Olp'^) fit, respectively. The different 
contributions to the corrections are also plotted, as well as the nonanalytic 
contributions -RAr,na (which are not included in Rn)- Note that the difference 
between the order p"^ and order p"^ curves stems only from (i^Q, which vanishes 
at order p"^. 

4.2 Systematic effects 

After having evaluated the size of the finite volume effects numerically, we 
briefly discuss effects which have been neglected in our analysis. There are 
three different types of systematic corrections: 

• We neglect the contributions from the nonanalytic part of the scattering 
amplitude, -RAr,na in Eq. (11). 

• Then there are the neglected higher order contributions of the Liischer 
formula which vanish as e"^*'^"^, with ^ > v3- 

• Furthermore, the pole of the A(1232) resonance is only taken into ac- 
count as a polynomial via the LEG in the subthreshold coefficients. 

In the following we argue why it is justified to completely neglect the first 
two of these effects and how we account for the last one. 

We show by an explicit calculation that the nonanalytic contributions 
are very small. It is possible to solve the dispersive integral Eq. (10) in 
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closed form, as described in Ref. [23]. Removing the fictitious poles from the 
expression 

Dfis) = ,,^Z .,3 {^t + al^t + 9\^t] /(^) (20) 



UF^tunv^ 



with 



\^) - 


1 / u^ f ^ \ 


SvrV Ml V Mj' 


^t - 


= 12z/^(-m^z/ + 3z/2-M2), 


At - 


= 2Au\u'-M'J, 


^t - 


= ^{v^-Mlf{-mMV + ^v^ + Ml 



leads to 



(21) 



Dt{s) = bt{s) - J2 dtn^'^- (22) 



The imaginary part of Df{s) vanishes identically in the Liischer integral, 
since it is an odd function in the integration variable. Furthermore, in for- 
ward scattering, the contribution from Df{u) is the same as the one from 
Df(s). Evaluating RN,na. numerically one finds that the contributions of 
the nonanalytic part of the scattering amplitude to the FVC can safely be 
neglected in the present analysis. 

Higher order terms in the Liischer formula for finite volume corrections 
to the pion mass have been explicitly calculated in Ref. [17]. It was found 
that for M^L > 2, the contributions which are not included in the resummed 
Liischer formula are very small. The formal order of the neglected corrections 
to the nucleon mass is the same as for the pion mass, e~^^^^. This suggests 
that the relative effects of these higher order terms for the nucleon mass are 
even smaller. 

Finally, to estimate the effect of the A(1232) resonance, the particle is 
included into chiral perturbation theory as an explicit degree of freedom. 
The resulting tree- level scattering amplitude -D^ can be found in Ref. [24]. 
Plugging this amplitude into the Liischer formula Eq. (2) we obtain R^. 

In Fig. 8, we compare R^ with Ra,cxp, which is the Liischer integral over 
the subthreshold expanded amplitude D^. The ratio P = (-Ra,cxp — Ra)/Ra 
is shown as a function of the quark mass for different box sizes. The following 
numerical values (from the set I fit, where relevant) are used for the nucleon 
and A(1232) masses and the A-nucleon coupling: 

rriN = m , m/^ = 1.232 GeV , qan = 13 . (23) 
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Figure 8: The plot shows the ratio P for different sizes of the box as a 
function of the pion mass. Note that although P grows fast for larger pion 
masses, the absolute correction R^ vanishes exponentially. The LEG from 
the order p'^ fit set I are used. In the region where the curves are dotted, 
M^L < 2. 



Here, a comment is in order. The additional pole of the A(1232) in the 
scattering amplitude might lead to modifications in the Liischer formula, as 
will be described in Sec. 5. For our choice of masses (see Eq. (23)) for the 
nucleon and the A(1232), this is however not the case. 

The corrections are large because for small values of M^^L, the suppression 
of the exponential factor in the Liischer integrand in Eq. (2) is not strong 
enough. The region with |z/|/M.^ > 1 — where the subthreshold expansion 
does not converge anymore — yields sizable contributions. The absolute size 
of Ra compares to a good approximation with the contribution of d^Q, but 
has the opposite sign. Therefore, the use of the subthreshold expansion in the 
Liischer formula leads to a systematic overestimation of the finite size effects 
if the pole of the A(1232) is not included explicitly in the vrA^ scattering 
amplitude. We account for this effect by introducing an asymmetric final 
error Rm ± at for the finite size effects. 
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The definition of a^^^ is provided in Appendix B. 



(24) 
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5 Resummed Liischer formula for the heavy 
meson masses 

In this section we present a formula that we will use for analyzing the finite 
volume corrections to the mass of a heavy meson similar to the one in Eq. (2). 
Again we use a resummed version of the Liischer formula. In the literature 
one finds a relativistic [30] as well as a nonrelativistic [31] effective Lagrangian 
that can be used to calculate the scattering amplitude. They are very briefly 
introduced in Appendix C, where we also give further references. 

Indeed for both formalisms we can provide a Liischer formula and then 
calculate the finite volume mass shift. We first discuss the formula that takes 
the relativistic amplitude as an input. The expression for the relative finite 
volume correction i?g' = {mH{L) — mH)/'fnH in this case reads as follows 



R 



rel 
H 



Ml 



E 



m(n] 



^ n=l ^ 



where 



M^ + ^j{ — "^H* 
2muM^ 



27re^glHH.e-V^^'-'')^- e{M^ + mj, - 

/oo 

-oo 

A, = M^L , f+{y) = T^iiM^y, 0). 



TTl 



H'j 



(25) 



(26) 



Here and later we adopt the notation that H and H* stand for any of the 
heavy pseudoscalar and vector mesons, respectively. The Liischer formula is 
valid provided that 



A* = tuh* — niH > 0, M^, A* ^ tuh, mn* 



K > 2. (27) 



The first two relations are well satisfied for all physical masses, while the 
third one requires L > 3 fm at the physical pion mass. T+ is part of the 
elastic ttH scattering amplitude, defined as 

T{7T\q)H{p) ^ 7r-\q')H{p')) = T,,, = S^'aT^ + \[ra'. rjT" , (28) 

where the amplitudes T^ depend on v and t, which are defined similarly as 
in Eq. (6): 



t= {q-q 



/\2 



S — U 



Am 



where s = {p + q)' 



u 



{p-q'y. (29) 



H 
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The proof of this formula is done in the same spirit as in the original 
article by Liischer [9]. The crucial difference is that the interaction of a 
heavy pseudoscalar meson H with a pion is mediated by the heavy vector 
meson H* , which has a slightly larger mass than the H . At low energy the 
pole generated by the if * exchange is the most important singularity and the 
one which dominates the Liischer formula. If one does not take the m^ — ?> oo 
limit and considers H* as nondegenerate with if, a third mass enters the 
FVC analysis and leads to some modifications of the proof by Liischer: the 
pole that is responsible for the first term in Eq. (25) is shifted outside of the 
integration contour for some configurations of the meson masses such that 
its contribution to the mass shift vanishes (see Fig. 1). This is accounted for 
by the step function. Note that the exponential in the first term can even 
become complex for some values of the masses, as the constraints given in 
Eq. (27) do not restrict e-^ to values smaller than L However, this only occurs 
for mass values where the step function vanishes such that the resulting mass 
shift remains real. 

We stress that in the degenerate case niH = rrifj*, Eq. (25) takes the 
same form (modulo trivial overall factors) as the one for the nucleon mass — 
which is the reason for discussing the two in the same paper. What makes 
this case particularly interesting is precisely the presence of a second almost 
degenerate mass, whose exact value we are free to vary. By varying it, we 
can better study and understand the role of the first, exponentially leading 
term in the Liischer formula, and ask, in particular, what happens if we 
change the masses such that the step function vanishes, and the first term in 
Eq. (25) disappears. This question is also relevant for the nucleon case if we 
consider the contribution of the A, which we have briefly discussed in Sec. 4.2. 
Indeed as we vary the pion mass, the A mass may get closer (depending on 
its quark mass dependence) to the nucleon mass and the corresponding pole 
may move from the right-hand to the left-hand side of the imaginary axis in 
the i/-plane. It is then a relevant question, what influence this has on the 
FVC. The dependence of the FVC on the mass difference A^, will be discussed 
below. 

The Liischer formula takes as input the amplitude in forward direction, 
i.e. for t = 0, and for u purely imaginary and small, as the contributions for 
large z/ are suppressed due to the exponential in the integrand. Similarly as 
for the vrA^ amplitude, it is convenient to use the subthreshold expansion of 
the ttH amplitude: 

T+(z/, 0) = Tp+ (z., 0) + T^{u, 0) + r+ (., u) , (30) 
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where 

,2 



^pl(->0) = ^f^ 



"^^' ' ^l-^'' (31) 

with 

M^ + mi, — mi,, , , 

vc = -£nM^ = '^-^ — ■ 32 

In contrast to the subthreshold expansion for the ttN amphtude, there is no 
inverse mass factor in the definition of T^, because here the amphtude is 
dimensionless. In order to have a dimensionless coupling constant Q-kHH* , we 
also need to apply a definition that is somewhat different from the usual one, 
namely 

2 2 

91hh* = lim^ '^^-j^ T . (33) 



Inserting the subthreshold expansion in the formula for the mass shift in 
Eq. (25), we get 



^ n=l * 

- t+ 5°(v^A^) + t+ 5'(v^A.) - it^B\^K)\ + RH,n. . (34) 

The finite volume integrals Jpv and B^ are the same as in Eq. (12). Further- 
more, i^Q = M^^t^Q and RH,na is the remainder coming from the analytically 
nontrivial part of the amplitude. 

At tree level, we find for the subthreshold coefficients 

.+ _ qIhh^^hi^c + 2M^)iyc .+ _.+ _ ^ .o.^ 

with the coupling constant given by 

Furthermore, Tj^(s, u) = 0, as no nonanalytical contributions can come from 
tree- level graphs, which then implies that -R^.na = 0. 
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The expression for the relative finite volume correction i?^'^' = {muiL) — 
mH)/mH in the nonrelativistic formalism reads 

RT = -TT^ f: ^ r dye-V^^^^^^f\y) . (37) 

Note that the absolute mass difference mnR^^^ is now independent of mu. 
The term proportional to q^hh* ^^^ disappeared because the argument of 
the step function is always negative if we take the limit tuh — ?■ oo at fixed 
A* (no matter how small it is). 

The proof of the Liischer formula only depends on the positions of the 
poles of the propagators, which are the same also in the nonrelativistic theory. 
Thus there is no need to repeat the full proof, all we have to do is to expand 
the absolute mass difference mnR'n'^^ in Eq. (25) to leading order in l/mn- 
As pointed out in Appendix C, the fields used for the two Lagrangians differ 
in normalization and to compensate for this we have to absorb a factor l/irtH 
into the scattering amplitude T"*". 

Because the nonrelativistic amplitude has mass dimension —1, we have 
again to change the definition of the coupling g-^HH* to 

A^ -zy2 
91hh* = lim, — ^T T , (38) 

which is dimensionless. This redefinition then leads to 



9ttHH*'^* 
A? -1.2 



rp>>0) = ^g^. (39) 



Inserting the subthreshold expansion in the formula for the finite volume 
effect in Eq. (37), we get 



pnrel _ ^l V" ^^^ 



-%^4V(v^A.,e.) 



327r2m// ^^ \/n\ 

n=l * 



RH,n.. (40) 



The finite volume integrals Jpv and B^ are again the same as in Eq. (12) and 

A, 



M. 



(41) 



Also, t^Q = M^HfQ and RH,na is the remainder coming from the analytically 
nontrivial part of the amplitude. 
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From the nonrelativistic Lagrangian we then get at tree level, 

91hh* = -4(A^ - M'^) , ^nt(^, «) = 0. (43) 

These results are identical to the leading term in the l/rriH expansion of the 
corresponding relativistic expressions (up to factors of niH and A* due to the 
normalization of the heavy meson field and the differing definition of q-khh* , 
respectively) . 



6 Numerical analysis of the FVC to the heavy 
meson masses 

The finite volume effects can now be calculated for arbitrary pion and heavy 
meson masses and box sizes L within the validity of chiral perturbation 
theory. To get meaningful results, one has to ensure that the conditions 
given in Eq. (27) are respected. 

For the coupling constant g a number of values from different sources are 
available in the literature (see Refs. [32-38]). g depends on both the mass of 
the heavy and the light quark, and neither dependence is fully understood. 
The listed publications give values for g for different heavy quark and light 
quark masses and, including also the error bars, these range from 0.18 to 0.79. 
For the numerical results presented in the following, we have used g = 0.5 
everywhere and thus neglected any quark mass dependence. 

Figures 9 and 10 show the finite volume effects for the D and B meson 
as a function of the pion mass for box sizes L = 2, 3 fm. The contributions 
coming from the Born term and from tJo are plotted separately for both, 
the relativistic and the nonrelativistic Liischer formula. The range of values 
covered by varying g from 0.18 to 0.79 in the relativistic formula is shaded 
in gray. Despite the large uncertainties the main message to be taken from 
the figures is that these FVC are negligibly small even for volumes as small 
as L = 2 fm. 

We find it nonetheless instructive to discuss in some detail some features 
which one can read off from the numerical analysis. First of all we remark that 
the deviation between the relativistic and the nonrelativistic result comes 
mainly from the t^Q term and it is very small compared to the uncertainty 
coming from the coupling constant. A detailed analysis of the size of the 
deviation will be given below. 
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Figure 9: Finite volume effects as a function of tlie pion mass for a box size 
of L = 2 fm, in the upper panel for the D, in the lower panel for the B meson. 
The contributions of the different terms to the Liischer formula are shown 
separately, the finite volume effect Rh is given by the sum. Solid lines are for 
the relativistic, dashed lines for the nonrelativistic formalism, the shaded area 
is the uncertainty on the relativistic result coming from the coupling constant 
g. Note that for L = 2 fm, the pion mass should fulfill M^^ > 0.2 GeV in order 
to stay in the region of validity of the Liischer formula. 



27 



0.0015 



0.001 



0.0005 



-0.0005 



0.0008 



0.0006 



0.0004 



0.0002 




M [GeV] 





1 




- 


— Rh 


- 




Bom 

-to; 


- 




^^--^ 






- 




1 







0.1 



0.2 



M [GeV] 



0.3 



0.4 



Figure 10: Finite volume effects as a function of the pion mass for a box 
size of L = 3fm: upper panel for the D, lower panel for the B meson. 
The contributions of the different terms to the Liischer formula are shown 
separately, the finite volume effect Rh is given by the sum. Solid lines are for 
the relativistic, dashed lines for the nonrelativistic formalism, the shaded area 
is the uncertainty on the relativistic result coming from the coupling constant 
g. Note that for L = 3fm, the pion mass should fulfill M^^ > 0.13 GeV in 
order to stay in the region of validity of the Liischer formula. 
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For the B meson, the Born term compensates parts of the very large 
t^Q contribution at small pion masses, while for larger pion masses, both 
contributions are positive and of similar size. For the D meson the Born 
term clearly dominates the finite volume effect. The net effect is smaller for 
the heavier B meson and the f^Q contribution is reduced more than the Born 
term, as it is suppressed by an additional factor of l/niH- 

In order to compare the results from the relativistic and the nonrelativistic 
framework we define the relative difference between the two as 

5h = R'^'/RT' - 1 • (44) 

Numerical evaluation shows that 6h indeed goes to zero for very large heavy 
meson masses. In the following, we give 6h for physical values of the masses 
for a box size of L = 3fm. For the D meson with rriD = 1.87 GeV and 
A, = 0.142 GeV it is about 0.75%, for the B meson with m^ = 5.28 GeV 
and A* = 0.046 GeV about 2.5%. Because Rh is proportional to g^, 6h is 
independent of the coupling constant and consequently we do not give an 
error on these numbers. The absolute value of 6h is very slowly decreasing 
with growing L such that the given values are good estimates for any box 
size. For unphysical pion masses, 6h can become quite large, indicating that 
the nonrelativistic framework is not very well suited for large pion masses. 
For Mtt = 0.4 GeV, say, it is about 10% for the B and about 17% for the D 
meson. 

Surprisingly, despite being the leading order term of a l/mn expansion, 
the nonrelativistic framework is more accurate for the lighter D meson. The 
reason is that 6h also depends on the mass splitting A*. If no heavy quark 
symmetry breaking effects are taken into account, i.e. when A^, = 0, the 
agreement is better for the B meson. 

We add some comments on the role of the step function in Eq. (34), which 
is nonzero if 



A, < A, = ^Jml + Ml-mH^ ^ . (45) 

By varying the parameter A,, between zero and the physical value observed 
for the case of the D and B mesons we can interpolate between the case in 
which the singularity structure of the amplitude is exactly like in the case 
of the nucleon (formally reached in the infinite quark mass limit), and the 
physical case with the H* meson nondengenerate with the H meson. This 
corresponds to moving the pole in the i/-plane from the left-hand side to the 
right-hand side of the imaginary axis, as shown in Fig. 1. How this affects 
the FVC is demonstrated in Fig. 11, where the finite volume effect is plotted 
as a function of the mass splitting A*. There one can clearly see the point 
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Figure 11: Finite volume effect for the D meson as a function of A^, witli 
L = 3 fm and M.^- = 0.4 GeV. Solid lines are for the relativistic, dashed 
lines for the nonrelativistic formalism. The Born term has been split up in 
the contribution containing the step function and the one with the integral 
Ipy. The latter has a discontinuity, which is compensated by the former that 
starts to contribute exactly at this point. The error band has been omitted. 

where the step function starts to contribute. At the same point, the other 
contribution to the Born term, containing the integral /pv, is discontinuous. 
This is due to the fact that e^r goes through zero at A^, = A^,, while the 
product ^TT-^pv goes to a nonzero value in the limit e,r — ^ and Ipy is even in 
En- The term proportional to the step function compensates this effect such 
that the resulting finite volume effect is continuous. Indeed this term arises 
from a deformation of the integration path which keeps the pole due to H* 
particle always to its right-hand side of the path, even when the pole moves 
to the left-hand side of the imaginary axis. Since in this way the integration 
path does not cross any singularity as we change A*, the integral remains 
analytic in A^, and we get a continuous curve. The figure also shows that in 
the nonrelativistic version of the Liischer formula the mass shift is continuous 
even without the step function term, because e.,^ is always negative for A^, > 
and thus the point of the sign change is not in the range of validity of the 
Liischer formula. 

Other than for the nucleon, where the physical observables entering the 
Liischer formula are best known at physical quark mass, the limited knowl- 
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edge on the coupling constant q^hh* does not allow for a better accuracy on 
the FVC at the physical point in the case of the heavy meson. Also, it is 
very challenging to perform lattice simulations for heavy mesons with very 
small masses for the light quarks and thus we refrain from discussing this 
situation in detail. 

7 Summary and conclusions 

In this paper we have evaluated the finite volume corrections to the mass of 
nucleons and heavy mesons in the p-regime on the basis of the resummed 
Liischer formula. The latter relates the finite size effect for the mass of a 
particle P to an integral over the forward irP scattering amplitude. We 
give simple analytic formulas which express the finite volume corrections to 
the nucleon and heavy meson masses in terms of only a handful of physical 
observables. 

In the case of the nucleon, these observables are the pion mass M,r, the 
proton mass niN, the pion-nucleon coupling constant g.„N and the subthresh- 
old parameters d^. Relying on the chiral representation of these quantities, 
one only needs to know the values of the chiral low-energy constants (LEG) 
to be able to predict the finite volume effects as a function of the quark mass. 
Fitting the nucleon mass rriN and the axial charge qa to lattice results as well 
as fixing the other observables at the physical point, we extract values for 
these LEG. These values are then used as a basis for our numerical analysis 
which shows that: 

• At the physical point, i.e. for a pion mass of M^. = 0.140 GeV, the 
relative finite size effects drop below 1% for L > 3.5 fm. 

• At the physical point and for box sizes of L > 3 fm we can predict the 
relative finite size effect with a precision of about l%o. 

• For unphysical pion masses, the relative finite size effects stay below 
10%. However, moving away from the physical point, the precision 
deteriorates quickly and an accurate prediction of the finite size effect 
becomes impossible. For instance, in a box of the size L = 2 fm and 
for Mjr = 0.3 GeV, the relative finite size effect can be between 0% and 
6%. 

• The finite size effects strongly depend on the order p'^ LEG Cj and on the 
quark mass dependence oi Qt^n (and therefore oi qa)- Future analyses of 
better lattice data should exploit these relations to extract interesting 
information about these quantities. 
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• Our final results for the relative FVC are tabulated in Table 2. 

We stress that it is thanks to the resummed Liischer formula that one can 
relatively easily evaluate higher order contributions to the final volume effects 
and better assess the uncertainties. At the plain one-loop level in chiral 
perturbation theory [11], for example, the quark mass dependence of qa does 
not yet play a role, and the contribution of the subthreshold coefficient d^Q 
is still absent. Both these effects are important, as we have shown, and 
substantially contribute to the uncertainties. 

In the case of a heavy meson much less is known about the scattering 
amplitude of pions off them. At present, the best we could do is to analyze 
this at tree level in the heavy meson chiral perturbation theory approach. At 
this order, the amplitude can be expressed in terms of the pion mass M^^ and 
the B^*^ and D^*^ meson masses as well as the coupling constant g describing 
the ttHH* coupling (with H = D or B). The numerical analysis is performed 
with these as input and shows that: 

• The limited knowledge on the coupling constant g leads to large error 
bars in the numerical results. A higher order calculation of the finite 
volume effect is only worthwhile once this has improved. 

• Despite our ignorance about the coupling constant g the FVC come out 
to be negligible (well below 1%) in all regions of the parameter space 
we have explored. 

Although the practical interest of a calculation of the FVC for heavy mesons 
appears to be limited (in view of the fact that they are small), we found a 
detailed analysis of some of its features very instructive. We have discussed, 
in particular, the difference between a relativistic and a nonrelativistic version 
of Liischer's formula and analyzed the dependence of the FVC on the size of 
the splitting between the pseudoscalar and the vector heavy meson. 
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M^[GeV] mAr[GeV] (5„^[GeV] Collaboration 



0.13957 


0.93827 


0.00008 


physical point 


0.192 


0.990 


0.037 


BMW 


0.208 


0.951 


0.021 


BMW 


0.271 


1.067 


0.030 


BMW 


0.307 


1.132 


0.031 


BMW 


0.318 


1.097 


0.026 


BMW 


0.320 


1.124 


0.026 


BMW 



Table 3: Data points for the fit of the nucleon mass, taken from Ref. [1]. 

A Low-energy constants 

The LEG c, and dis are scale independent. For the remaining LEG the scale 
independent combinations are given by 

jS p ^m fi 

^, , 1 + 3^2 + f ^^ + 8cim + Cam - 4c3m , M^ 
es = e^ /i + ^ -— In , 

10 + 12g^ + fg^ + 8c2m M, 
'' - '^^^^ NF^ ^^~' 

ee = e,{^)+ ^^,^ In—, (46) 

with N = 1671^. 

B Details of the fits 

The lattice results for the nucleon mass and the axial charge which are used 
in the fit are shown in Tables 3 and 4. We intend not to strain the chiral 
expansions too much and do not include lattice results with pion masses 
significantly above 350 MeV. For the remaining quantities, only the value 
at the physical point, known from experiment [42], is available, 

g^N = 13.39 ± 0.08 , M^doo = -1-46 ± 0.10 , 
M^dio = 1.14 ± 0.02 , M^d2o = 0.200 ± 0.005 . (47) 
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5.. 


Collaboration 


0.13957 


1.269 


0.003 


physical point 


0.313 


1.230 


0.100 


ETM 


0.350 


1.210 


0.070 


RBC/UKQCD 


0.352 


1.250 


0.060 


LHPC 



Table 4: Data points for the fit of the axial charge. The data points are 
taken from Refs. [39-41] 



The LEG ii from the mesonic order p^ Lagrangian are not determined in the 
fit, but are taken as input parameters without an error. We use the values 
given in Ref. [43]. For the pion decay constant F,^, the two-loop result from 
Ref. [44] with the values of the LEG as discussed in Ref. [10] is used. We use 
the standard x^ function 

i,k ^ ' ^ 

where yk{'x,Mi) is the chiral representation of quantity i/k evaluated at the 
pion mass Mj. The LEG are collected in the vector x and yk^i denotes the 
data point of quantity yk at the pion mass Mj with the error 6k^i- The error 
(Tyj. of the quantity yk is given by 



with C the correlation matrix of the LEG. 

C Formalism of heavy meson ChPT 

At present there are two frameworks to include heavy mesons into chiral 
perturbation theory: a relativistic one introduced in Ref. [30] and a nonrel- 
ativistic one introduced in Ref. [31]. We only present the two Lagrangians 
very briefly; for a detailed discussion, we refer to the original publications 
and to Refs. [45, 46]. All conventions we use have been adapted to agree 
with the latter. We give the relativistic Lagrangian for fields Ba and Ba 
that annihilate a pseudoscalar and a vector meson consisting, respectively, 
of a heavy quark Q and a light antiquark g^. From these we build the field 



lo I i^ab + f^H^ab] r . j-, , j-,*/i i ,0 Oct 

" " I 2^ ) ^ '''^^ ^b^fA^ ^^ = T ^I 



l7°, (50) 
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where we omitted the superscript for simphcity. The covariant derivative is 
given by D^^^ = d^5ab — T^^, with the connection 

r^ = ]^[u\ d,u] = \{u^d,u + ud^u^) = -^ e'^^W'a^ + 0{tx^). (51) 

We also need the vielbein 

u^ = i{u\ d^u} = i{u^d^u - ud^u^) = -4^"9^vr» + (9(7r3). (52) 

r 

From these building blocks we can construct the relativistic heavy meson 
Lagrangian, which is 



C-HM, rel — —rRH Tr 



Ba{0ab - mHSab)Bb - -BaBbJ^.J^u'' • (53) 



The second term is not contained in the Lagrangian of Ref. [30] and had to 
be added here as it is the source of the one-loop self-energy contributions. 

By means of a nonrelativistic reduction procedure similar to the one used 
in heavy-baryon chiral perturbation theory [47, 48], we can derive the La- 
grangian of Ref. [31] from this. Instead of Ba, we use the field 

Ha=^^ [^Pal, + P:^%] , H^ = I'Hl^' = [zPi75 + P7^7.] ^ ^ (54) 

where f ^ is the four velocity of the heavy meson and the fields Pa and P*'^ 
replace Ba and B^^ respectively. Note that the new fields differ in normaliza- 
tion by a factor l/y/rriH, which ensures that the nonrelativistic Lagrangian 
does not depend on rriH- Inserting these fields into Eq. (53) and expanding 
to leading order in l/ruH, we find the nonrelativistic Lagrangian 



ChM, nrel = -i Tr [HaV^d^Ha) 
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+ ^ Tr {Ha H,v,Tl) + | Tr {5^ H,^,^Xa) ■ (55) 

Up to the phase conventions, this agrees with the nonrelativistic Lagrangian 
of Ref. [31]. 
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